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ABSTRACT 

We study the high-energy emission of the Galactic black hole candidate GX 339^ using INTEGRAL/SPI and 
simultaneous RXTEfPCA data. By the end of January 2007, when it reached its peak luminosity in hard X-rays, 
the source was in a bright hard state. The SPI data from this period show a good signal to noise ratio, allowing 
a detailed study of the spectral energy distribution up to several hundred ke V. As a main result, we report on the 
detection of a variable hard spectral feature (>150keV) which represents a significant excess with respect to 
the cutoff power law shape of the spectrum. The SPI data suggest that the intensity of this feature is positively 
correlated with the 25 - 50 ke V luminosity of the source and the associated variability time scale is shorter than 
7 hours. The simultaneous PCA data, however, show no significant change in the spectral shape, indicating that 
the source is not undergoing a canonical state transition. We analyzed the broad band spectra in the lights of 
several physical models, assuming different heating mechanisms and properties of the Comptonizing plasma. 
For the first time, we performed quantitative model fitting with the new versatile Comptonization code belm, 
accounting self-consistently for the presence of a magnetic field. We show that a magnetized medium subject 
to pure non-thermal electron acceleration provides a framework for a physically consistent interpretation of 
the observed 4 - 500 keV emission. Moreover, we find that the spectral variability might be triggered by the 
variations of only one physical parameter, namely the magnetic field strength. Therefore, it appears that the 
magnetic field is likely to be a key parameter in the production of the Comptonized hard X-ray emission. 
Subject headings: methods: observational - X-rays: individual (GX 339^) - accretion, accretion disks - 
radiation mechanisms: general - magnetic fields 



1. INTRODUCTION 

GX 339-4 was discovered in the earl y 70's by the MIT X - 
ray detector aboard the OSO-7 mission dMarkert et alj|1973b . 
The source is classifi ed as a low mass X-ray binary (LMXB; 
IShahbaz et alj|200l1) from upper limits on the optical lumi- 
nosity of the c ompanion star. It is believed to harbor a black 
hole, for which lHynes et alJ (l2003h derived a mass function of 
5.8 + 0.5 Mq. The inclination o f the system is yet unc ertain. A 
study of the binary parameters ( IZdziarski et alj2004h revealed 
a plausible lower limit of / > 45°. while ICowlev et al.l ( l2002h 
suggested / < 60° because the system is not eclipsing. On 
the other hand, spectral fits of the Fe K„ region with Chandra 
JMiller et al. 2004b) and XMM Newton (Miller et al. 2004a; 
iReis et ani2008i) clearly favor lower values for the inner disc 
inclination (typically / ~ 20°). This may indicate that the 
inner accretion disc is warped, making it appear at a lower in- 
clination than the orbital plane. Similarly, there is no certainty 
concerning the distance to the source. Reso l ving t he velocity 
structure along the line of sight. iHynes et ar(2004) obtained a 
conservative lower limit of d > 6kpc. IZdziarski et al. (2004) 
analyzed the binary parameters of the system and found the 
most plausible distance to be around Skpc. In this work, we 
use d - Skpc, / = 50° and, assuming a small companion 
mass, we infer M = 13 Mq from the mass function. 

Black hole binaries (BHBs) are powerful engines produc- 
ing high-energy radiation up to the y-ray domain. From 
an observational perspective, they appear in four spec- 
tral state s, namely, quiescent, low/hard, intermed iate and 
high/soft ("Tanaka & Lewin"1995L lBelloni et al.ll2005l: see also 
iMcCIintock & Remillard.200a for a slightly different classifi- 
cation). 



Soft states: In the soft state, the energy spectrum is domi- 
nated by a soft (< lOkeV) component which is attributed to 
thermal emissi on from an optically thick geometrically thin 
accretion disc (IShakura & Sunyaevlll973h extending down to 
the last stable orbit. Above lOkeV, the emission is charac- 
terized by a complex hard X-ray continuu m with no clear 
presenc e of a high- e nergy cutoff (see e.g. iGierliriski et alj 
(1999); iMotta et al.1 (l2009h : ICaballero-Garcia et al.l (l20'09la 
This component can be attributed to inverse Compton scat- 
tering of soft photons (UV, soft X) in a hybrid thermal/non- 
thermal electro n plasma (the so-called 'corona') (see e.g. 
the review by iDone. Gierlihski. & Kubotal 120071 hereafter 
DGK07). 

The non-thermal particles in the corona may be gener- 
ated due to the magnetic field. Indeed, the Parker insta- 
bility is able to transport a significant fraction of the ac- 
cretion power above and belo w the disc (iGaleev et al.l 1 19791 ; 
lUzdenskv & GoodmanI |2008|) . where the energy may then 
be dissipated through magnetic reconnection. Alternatively, 
Fermi acceleration at relativistic shocks is also expected to 
produce non-thermal particle distributions which can explain 
the observed steep power law emission in hard X-rays. 

Hard states: In the hard state, the energy spectrum is very 
different. The disc emission almost vanishes and the spec- 
trum is dominated by a hard power law component (photon 
index F - 1.4-2.0) with a nearly exponential cutoff (Zscut - 
50 - 150keV) (see e.g. [Zdziarski et al. 1998). Such a spec- 
trum is well described by thermal Compto nization in a hot, 
optically thin electron-proton plasm a (iSunyaev & Truemped 
1 1 9791: IZdziarski & Gierlihskill2004l) . Soft y-ray observations 
of several BHBs additionally revealed a high-energy excess 
with respect to a thermal Comtonization model, suggesting 



that in some cases at least some level of non-thermal electron 
acceleration is also required. For instance, COMPTEL ob- 
served a non-th ermal tail in the avera ged hard state spectrum 
of Cygnus X-1 dMcConnell et al.l |2002) while OSSE and SPI 
de tected such a feature during bright hard states of G X 339- 
4 dJohnson et all 119931: IWai-dziriski et al. 2002; Joinet et alj 
I2OO7I) ! 

Intermediate states: Intermediate states are observed dur- 
ing transitions between the hard and the sof t states and usu- 
ally show the characteristic features of both ('Mivamot o et alj 
[1991; Men dez & van der Klis 19 97; Kong et al. 2002). For 
GX 339-4, iBelloni et aH ( 120051) defined 'two diff'erent vari- 
eties of intermediate state: the hard intermediate state (HIMS) 
and the soft intermediate state (SIMS). The transition from the 
HIMS to the SIMS of th e 2004 outburst of GX 339-4 is re- 
ported in lDel Santo et al.l (120081) . 

According to a popular scenario, the different spectral states 
can be explained through changes in the geometry of the ac- 
cretion flow. The weakness of the thermal component in the 
hard state is generally interpreted as a consequence of a trun- 
cated accretion disc (DGK07), which, in its inner parts, is re- 
placed by a hot, advection dominated accretion flow (ADAF; 
IShapiro et alTTgv e; Narayan & Yi 1994; Yuan et al. 2007). In 
these solutions, gravitational energy is converted into thermal 
energy of protons, which in turn heat the electrons through 
Coulomb collisions. This process naturally forms the quasi 
thermal electron distributions that are required to explain the 
typical hard state spectra. Moreover, recent models of hot ac- 
cretion flows include a small non-thermal component that can 
account for the sometimes observed high-energy excess. 

However, a number of recent results seem to question 
this paradigm. First, it is notoriously difficult to estimate 
the inner disc radius in the har d state and its recess ion is 
a highly debated topic (see e.g. ICabanac et al.]l2009h . Al- 
though several reports support the truncated disc scenario 
(iTomsick et al.ll2009l: lOone & Diaz Trigoll2009l) . others sug- 
gest that even in hard states the accret ion disc may ex- 
tend down to a few Schwarzschild radii dMiller et al.l 120061 ; 
iReis et ani2008h , which would be inconsistent with the stan- 
dard ADAF scenario. Second, hot accretion flows generally 
exist only if the optical depth is small {tj ^1), which im- 
plies that the electron-proton coupling is weak and therefore 
allows high proton temperatures. However, it is difficult to 
align such a small optical depth with both the hard X-ray 
spectral slope and the thermal cutoff energy. When the spec- 
tral slope is reproduced, classic ADAF models lead to higher 
electron temperatures than those obs erved in BHB s, as it was 
shown for Cygnus X-1 or GX 339-4 (lYuan & Zdzi arski 200l 
iMalzac & Belmont! l2009l hereafter MB09). Recent Monte- 
Carlo simulations, which account for global Com pton cool- 
ing and general relativistic effects dXie et al.l20IOl) . yield self- 
consistent solutions in which the peak energy of the spectrum 
is reduced with respect to previous results. However, it is still 
not straightforward to accommodate the inferred high-energy 
cutoff to the observations, even if gravitational redshift is ac- 
counted for 

As an alternative to the ADAF-like models, the Comptoniz- 
ing medium in the hard state could as well be powered by 
the same non-thermal mechanisms that are believed to ac- 
celerate the electrons in the soft state, i.e. diffusive shock 
acceleration or magnetic reconnection. Such a model natu- 
rally accounts for the presence of a non-thermal component 
in the high-energy spectrum. In addition, MB09 showed 
that the steady state electron distribution can appear quasi 
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Figure 1. Overview of the 2007 outburst showing the daily light curves of 
GX 339^ obtained by SWIFT/BAT (1 5 - 50 keV ; top panel) and RXTE/ASM 
(1.3-12 keV ; lower panel). The vertical dotted lines indicate the time period 
during which the here presented SPI data were obtained. 

thermal even if acceleration mechanisms are purely non- 
thermal (see also iPoutanen & VurnJ ( 120091) ). These authors 
studied the thermalizi ng effects of the magnetic field , since 
it was pointed out by iGhisellini. Haardt. & SvenssonI (Il998h 
that the very fast emission and absorption of synchrotron 
photons (the so-called 'synchrotron boiler effect') is able to 
thermalize the electron distribution in a few light-crossing 
times. Using belm, a new code whic h includes this effect 
dBelmont. Malzac. & Marc owith 2008!), MB09 qualitatively 
explained the variety of spectral states observed in the proto- 
typical BHB Cygnus X-1. The model is consistent w ith a disc 
recession in the hard state (IPoutanen & Vurnj2009l) . but since 
the weakness of the soft co mponent could as w e ll result from 
a low er disc temperature (iBeloborodovl 1 1 999t IMalzac et al.l 
12001b . it does not necessarily require a change in the geome- 
try of the accretion flow (MB09). 

In this paper, we use observations of GX 339-4 with 
SPl/INTEGRAL and PCA/RXTE to study the spectral energy 
distribution of the bright hard state from the 2007 outburst. 
We analyze the spectral variability in the framework of the 
models cited herebefore. The data and the associated reduc- 
tion procedures are described in Section 2. Section 3 is ded- 
icated to a phenomenological analysis of the high-energy be- 
havior while in Section 4 we present an extensive physical 
analysis of the broad band spectra. The results and their im- 
phcations are discussed in Section 5 and summarized in the 
conclusive section of the paper 

2. OBSERVATIONS 

The here presented data on GX 339-4 were obtained with 
the INTEGRAL and RXTE observatories at the end of January 
2007. At that time, the source was ver y bright in hard X-rays 
(cf. Figure [1) and iMotta et akl (120091) identified its spectro- 
temporal properties to be characteristic of the hard state. 

2.1. Instruments 

Our analysi s focuses on th e results obtained by the spec- 
trometer SPI (iVedrenne et al.l l2003), which is one of the two 
main instruments aboard the international gamma ray obser- 
vatory INTEGRAL. Operating in the 20keV - 8MeV band, 
SPI uses germanium detectors to provide high spectral res- 
olution while imaging is performed through the coded mask 
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Instruments 


Obs. ID 


MJD start 


MJD stop 


Exp 


time (ks) 


SPI/IBIS 


525 


54130.6 


54132.2 




107 


SPI/IBIS 


low (cutoft) 


— 


— 




36 


SPI/IBIS 


high (excess) 


— 


— 




36 


PCA/HEXTE 


92035-01-01-02 


54131.10 


54131.17 




3.7 


PCA/HEXTE 


92035-01-01-04 


54132.08 


54132.15 




3.7 



Table 1 

Log of the different observations discussed in the paper. 

technique (Roques et al. 2003). The observational strategy 
of the INTEGRAL mission is to sample the region of inter- 
est by means of 30 - 40 minute long fixed pointings (the so- 
called sci ence-window s ), each separated by a 2° angular dis- 
tance (see lJensen et alj (I2003h for details). In order to extend 
the spectral coverage to lower energies, we also considered 
simultaneous 4-28 ke V data from the P roportional Counter 
Array (PCA) on RXTE jBradt et al.lll"993h . 

2.2. Data and reduction methods 

The SPI data from the GX 339-4 region were obtained dur- 
ing INTEGRAL revolution 525, lasting from 2007 January 
30 to February 1 (cf. Table 1). We selected the science- 
windows where GX 339-4 was less than 12° off the cen- 
tral axis and which did not show any contamination by so- 
lar flares or radiation belt exit/entry. This resulted in a set 
of 50 useful pointings representing 107ks of observational 
coverage. In order to determine the emitting sources in the 
field of view, w e performed 25 - 50keV imaging with the 
spiRos software (ISkinner & Connenil2003h . Aside from the 
main source, 4U 1700-377, OAO 1657-415, IGR J16318- 
4848 and GX 340+0 were detected above a 5 cr threshold. The 
positions of these sources were then given as a priori infor- 
mation to a specific flux-extraction algorithm, using the SPI 
instrument response for sky-model fitting. In a first run of the 
software, we allowed the background normalization as well 
as the source fluxes to vary between successive pointings. In 
this way we determined the most appropriate variability time 
scale for each component. The background normalization was 
more or less stable during the SPI observation, therefore we 
assumed no background variability in the final reduction pro- 
cess. Some of the sources in the field, however, showed con- 
siderable flux variations. Most notably, 4U 1700-377 was ex- 
tremely variable, with a 20 - 50 keV flux per science-window 
ranging from to 700 mCrab. 

There are two RXTE observations which are simultaneous 
with the INTEGRAL exposures (see Table 1). The first ob- 
servation took place towards the middle and the second to- 
wards the end of revolution 525 (cf. the shaded areas in Fig- 
ure |2]i. For both observations, we downloaded the standard 
products from the HEASARC w ebsiteQ and analyzed the data 
in xsPEC vll.3.2 (lArnaudlll996b . Apart from a global 3 per 
cent flux difference, the PCA spectra from both observations 
are compatible within the error bars and were co-added us- 
ing the addspec ro utine from the ftools package. Following 
iMottaet aH (l2009h . we added a systematic error of 0.6 per 
cent to each PCA channel. 

2.3. SPI light-curve and data subsets 

Figure |2] shows the 25 - 50keV SPI light curve from GX 
339-4 obtained during revolution 525. Each bin represents 

' http://heasarc.gsfc.nasa.gov/docs/archive.htmI 
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Figure 2. 25 - 50keV SPI light curve from revolution 525. The pointings 
selected for the high and low flux data sets are shown in red and blue, respec- 
tively. The cyan line represents the average flux from the SPI observation. 
The shaded areas indicate the time periods of the simultaneous PCA obser- 
vations. 

one pointing or science-window, which corresponds to a time 
scale of about 40 min. Overall, the 25 - 50keV flux shows 
no particular evolution, indicating that the source remained 
in the bright hard state during the ~1.5 day long observation. 
However, the light-curve reveals minor variability (on a time 
scale of hours, or less) around the average flux value {<F>- 
657.1 + 3.8mCrab). To investigate whether this variability 
could be linked with changes in the spectral energy distribu- 
tion, we fixed two arbitrary flux bounds (Flow - 640 mCrab 
and Fhigh = 670mCrab) and grouped the science-windows 
where Fscw < ^low and Fscw > ^high, respectively. For each of 
the two resulting data subsets, color-coded in blue and red and 
referred to as low and high respectively, we then produced the 
averaged 25 - 500 ke V spectrum. Both subsets consist of an 
equal number of science-windows, which allows a straight- 
forward comparison of the spectra. Since the source flux is 
outstandingly high, the spectral extraction yields significant 
results, even if the number of science-windows is relatively 
small. 

3. HIGH ENERGY EMISSION 

3.1. SPI spectra 

The resulting 25 - 500 keV spectra differ by ~14 per cent 
in terms of their total flux and were separately fitted in xspec 
with a simple cutoff power law model. The uncertainty on the 
model parameters is given at the 90 per cent confidence level 
{l^x^ - 2.7). As the SPI data alone do not allow to simulta- 
neously constrain the photon index and the cutoff energy, we 
fixed the former to F = 1.3 and left the latter free to vary. The 
model provides a good fit to the low spectrum {x^ 122 - 0.75) 
and we infer a significant cutoff at E^ - 58.6 + 2.2 keV. 
For the high spectrum, the marked curvature around 40 ke V 
still suggests the presence of a cutoff energy, fitted at E^ - 
56.2 ±2.1 keV. With respect to the model, however, we ob- 
serve a significant excess above 150 keV, leading to an unac- 
ceptable quality of the fit Ot^/28 = 1.82). To account for the 
high-energy excess, we added a second power law of fixed 
photon index F = 1 .(fl The presence of this second compo- 

^ If this parameter is left free to vary, the fitted 90 per cent confidence 
interval is re [0.3,1.9] 
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Figures. Top panel: 150^50 keV SPI lightcurve with a binning of 7 h. The 
flux is given in units of 10^ counts s"'. The cyan cross indicates the average 
flux and flux uncertainty. Middle panel: [150^50]/[50-150] keV hai'dness 
ratio as a function of time. Lower panel: [150-450]/[50-150] keV hardness 
ratio as a function of the 25 - 50 keV flux. 

nent, which reduces the cutoff energy to Ec = 49.5 + 3.8 keV 
but improves the fit to y^/21 = 1.03, is statistically required 
as shown by the ftest (iBevington & Robinsonll992h . Indeed, 
we infer a probability of Pftest = 6.1 x 10"^ that the im- 
provement of the fit was a chance event. As a consequence, 
we conclude that both SPI spectra essentially differ in terms 
of the highest energy emission (> 150 keV). 

3.2. Imaging and spectral robustness 

As the flux extraction process of the SPI data can be sensi- 
tive to background modelling and a priori information about 
the distribution of emitting sources in the field of view, we 
double checked our results. In the 150-450keV spiros image 
drawn from the high data set, GX 339-4 is detected without 
any a priori information at a flux level of 347.4 + 49.9 mCrab 
and a 7.0 cr significance. This contrasts to a non-detection in 
the low data set, where the flux significance at the source po- 
sition is below 2.5 a. We note that the maximum level of the 
uniformly distributed residuals is around 4 cr in both images. 

Since GX 339^ is the only source detected above 150keV, 
we used a sky model with a single point source and re- 
extracted the high-energy part (>150 keV) of the spectra. We 
also tested different background models along with various 
variability time scales. The obtained results are all perfectly 
compatible within the 1 cr errors, showing that the spectra are 
insensitive to the details of the flux extraction process. We 
conclude that the respective presence of a cutoff and a high- 
energy excess in the SPI data is robust and significant. As a 
consequence, we rename the low and high spectra to cutoff 
and excess, respectively. 

3.3. Time evolution 

The presence and absence of a high-energy excess in the 
high and low flux spectra suggests that the appearance of this 
feature is correlated with the 25 - 50 keV flux of the source. 
To inverstigate this issue further, we extracted light curves in 
the 25 - 50, 50-150 and 150 - 450 keV energy bands with a 



binning of 10 science-windows (»7h) and plotted the [150 - 
450]/[50- 150] keV hardness ratio as a function of time and as 
a function of the 25 - 50 keV flux. From these plots, shown in 
Figure[3j it is clear that the strength of the high-energy excess 
is correlated with the 25 - 50 keV flux. Since a different time 
binning (8 or 12 science-windows for instance) leads to the 
same results and since the 25 - 50 ke V flux is a good tracer of 
the total X-ray luminosity, we conclude that the intensity of 
the high-energy tail is most likely correlated to the total X-ray 
luminosity of the source. 

For statistical resons, the SPI data do not allow to quantify 
the exact time scale of the phenomenon. Nevertheless, from 
the 25 - 50keV science- window light curve (cf. Figure |2]i 
and the above analysis, we can assess that the time scale of 
the observed spectral variability is shorter than 7 h. 

3.4. Cross-check with other instruments 



Instrument 


Model 


ar'/v)i 


(x^/y)i 


PpTEST 


SPI 


CUTOFFPL 
EQPAIR 


36/26 
47/26 


18/25 
22/25 


4x10-5 
2 X 10-5 


IBIS/ISGRI 


CUTOFFPL 
EQPAIR 


32/29 
39/29 


26/28 
25/28 


2x10-2 
5 X 10-* 


HEXTE 


CUTOFFPL 
EQPAIR 


59/44 
66/44 


47/43 
47/43 


2 X 10-3 
1 X 10-* 


SPI+IBIS+HEXTE 


CUTOFFPL 
EQPAIR 


145/103 
141/103 


105/102 
99/102 


1 X 10-** 

2 X 10-' 



Table 2 

Quantification of the high-energy excess in the total averaged SPI, 
IBIS/ISGRI and HEXTE spectra with respect to a cutoff power law and a 
thermal Comptonization model. To account for the high-energy excess, we 
added a second power law component with a fixed spectral index of F = 2.0 
in the phenomenological model, while we allowed for non-thermal heating 
in the physical model. Each time we give the fit quahty before (x^ /v)i and 
after (x'/v){ adding a degree of freedom, along with the ftest probability 
that the fit improvement was a chance event. 



Since the detection of a high-energy tail is a critical is- 
sue, we cross-checked the SPI results with other instruments. 
We analyzed the data obtained simul taneously by the soft 
gamma ray imager IBIS/ISGRI (lUbertini et al.ll2003l) aboard 
INTEGRAL and the high energy X-ray timing experiment 
HEXTE aboard RXTE. The HEXTE spectra from the two si- 
multaneo us RXTE observatio ns (cf. Table 1) have been pre- 
sented bv lMottaet al.l (120091) . who kindly provided us the re- 
duced data. Both spectra are compatible within the Icr er- 
rors and were co-added to obtain better statistics at high en- 
ergies. The total SPI and IBIS/ISGRI spectra (averaged over 
the whole INTEGRAL revolution 525) have been presented 
by ICaballero-Garcia et al.l (120091) . who report on the detec- 
tion of a high-energy excess above 150 keV. We re-extracted 
the IBIS/ISGRI data using the standard OS A 8.0 software 
packageQ and jointly fitted the averaged SPI, IBIS/ISGRI and 
HEXTE spectra in xspec. We find a good agreement between 
the three instruments, not only in spectral shape but as well in 
normalization. In addition, the data from all three instruments 
confirm the presence of a high-energy excess with respect to 
a phenomenological cutoff power law and a thermal Comp- 
tonization model (cf. Figure|4] Table 2). For the latter we used 

3 available from the INTEGRAL Science Data Centre (ISDC) 
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EQPAiR ( ICoppilll999l cf. section 4.2) with non-thermal heat- 
ing turned off (/mh = 0). The reflection amplitude is poorly 
constrained and was therefore fixed to Q/2:7r = 0.35 (average 
value derived from the physical analysis, cf. Table 3). 

The two IBIS/ISGRI spectra averaged separately over the 
low and high data sets are consistent with the SPI results. 
However, due to the shorter exposure times, the high-energy 
excess is no longer required in the ISGRI data. This is not sur- 
prising since the sensitivity of ISGRI above 200 ke V is lower 
than the sensitivity of SPI. As a consequence, we did not con- 
sider the IBIS/ISGRI data in the remainder of the paper be- 
cause they do not improve the high-energy constraints with 
respect to the models. The HEXTE data, on the other hand, 
are obtained with even shorter exposure times and are thus 
unsuitable for our group-wise spectral analysis. 

4. BROADBAND SPECTRAL ANALYSIS 

In order to understand the processes that are likely to cause 
the observed X/y-ray emission, we analyzed the broad band 
PCA/SPI spectra in the lights of different physical Comp- 
tonization models. Since the spectral variability is occurring 
at high energies only, each one of the two SPI spectra (cf. 
section 3.1) was jointly fitted with the co-added PCA spec- 
trum. A variable multiplicative factor was added to account 
for cross-calibration uncertainties as well as the luminosity 
difference between the two SPI spectra. The difference to 
unity of this factor is never larger than 8 per cent. Except for 
slight changes in the x^ values, neither the qualitative nor the 
quantitative results are affected when using one of the single 
PCA spectra instead of the co-added spectrum. 

4.1. General model 

The high-energy source in GX 339-4 is modeled by a 
spherical, magnetized, fully ionized proton-electron/positron 
plasma of radius R. The emission of the plasma is de- 
rived by self-consistent computations of the equilibrium elec- 
tron distribution accounting for Compton scattering, syn- 
chrotron emission/absorption, pair production/annihilation 
and Coulomb collisions {e-e and e-p). The Thomson optical 
depth of the plasma is given by tj = Tion + Tp^h-, where rion 
represents the optical depth of ionization electrons and Tpai, is 
the opacity that arises from pair production. 

The plasma properties essentially depend on the magnetic 
field strength as well as on the power supplied by external 
sources. The tangled magnetic field strength is parameterized 
by the magnetic compactness: 



Ib 






(1) 



Energy injection is quantified by the compactness parameter: 



/ = 



a-j L 



nieC^ R ' 



(2) 



where L is the power supplied to the plasma, cr-j- the Thom- 
son cross-section and mg the electron rest mass. The general 
model comprises three possible channels for providing energy 
to the coupled electron-photon system: (i) non-thermal elec- 
tron acceleration (/nth), (ii) thermal heating of the electron dis- 
tribution (/(h) and (iii) external soft radiation from a geomet- 
rically thin accretion disc (/s). The non-thermal acceleration 
processes are mimicked by continuous electron injection at a 
power law rate with index Finj (i.e. «(y) oc y^i^ini) and Lorentz 
factors ranging from ymin to ymax- The thermal heating of the 



electron distribution could be caused by Coulomb interactions 
with a population of hot ions. The incident radiation from 
an accretion disc is modeled by a blackbody component of 
temperature kTbi,. Since the temperature of this component 
cannot be constrained by our data, we follow iDel Santo et al.l 
(i2008b and fix kT^h to the fiducial value of 300 eV. All the 
injected energy ends up into radiation, and in the compact- 
ness formalism the total radiated power in steady state holds: 

I = 'nth + 'th + h- 

For a given source flux F, the total compactness can be es- 
timated by the formula: 

where M is the mass of the black hole and d the distance to the 
source. Here we expressed / in terms of Fq - 2.6 x 10"^ erg 
s^'cm"^, the observed average 4 - 500 keV flux of GX 339- 
4. Assuming a spherical plasma of radius R = 30 Rg, a black 
hole mass of M = 13Mo and a distance of c/ = Skpc, this 
leads to a total compactness of / ~ 100. 

Besides being a source of soft seed photons, the accretion 
disc may also give rise to Compton reflection of the inci- 
dent hard X-rays from the Comptonizing plasma. The re- 
flected emission is calc ulated using the viewing-angl e depen- 
dent Green's functions (Ma gdziarz & Zdziarskil 19951) . but ac- 
counting for general relativistic effects. The spatial extension 
of the possibly ionized disc is parameterized by its inner and 
outer radii, fixed at Ri„ = 6 Rg and Rout - 400 Rg- 

Moreover, the model takes into accou nt the Kg fluorescen ce 
of the Fe elements in the disc (diskline, iFabian et al.llI989() as 
well a s the inters t ellar a bsorption (phabs). For the latter, we 
follow lReis et alJ (l2008h and assume a fixed neutral hydrogen 
column density of A^h - 5.2 x 10"' cm"^. 

Since our data sets are energy and sensitivity limited and 
since such a complex modelling is partly degenerate, spectral 
fitting is unable to provide simultaneous constraints to all the 
parameters. Therefore, starting from the general model, we 
define several canonical sub-models where some of the pa- 
rameters are held invariant. 

4.2. Thermal heating 

Given that the relative ratio of the heating mechanisms 
(thermal versus non-thermal) is very difficult to constrain, we 
investigated the extreme situations in which there is only one 
possible channel for providing energy to the electrons. 

First, we consider the case where the energy injection is 
purely thermal (/th ^ 0, /nth = 0). The seed photons may 
arise either from an external source (i.e. the accretion disc), 
or, in presence of a magnetic field, from internally generated 
synchrotron emission. However, since the magnetic field is 
believed to give rise to non-thermal electron acceleration, we 
limited our analysis to the standard case where magnetic pro- 
cesses are neglected (Ib - 0). 

The main model parameters are the power supplied to the 
Comptonizing electrons /th, the power contained in the soft 
photons from the disc /,, the Thomson optical depth of the 
plasma rion, the reflection amplitude H/ln, the ionization 
parameter of the reflecting material ^ and the normalization 
factor between the two instruments. To calculate the model 
spectra, we use th e hybrid therm al/non-thermal Comptoniza- 
tion code eqpair JCoppil [19991) . A detailed description of 
the code and its application to Cygnus X-1 can be found in 
,Gierliriskietal..a999,) . 





Figure 4. Total averaged SPI (red), IBIS/ISGRI (blue) and HEXTE (green) spectra, jointly fitted with a thermal (left) and a hybrid thermal/non-thermal (light) 
Comptonization model. The normalization was fixed to the calibration of SPI and we applied mutiplicative constants of Cisgri = 1-03 and Chexte = 1 ■ 1 to the 
IBIS/ISGRI and HEXTE spectra, respectively. The lower panels indicate the data/model ratios for each instrument separately. These show that a pure thermal 
Comptonization model is not able to account for the emission above 150keV. 



Since the model spectrum strongly depends on the ratio 
khih, but only weakly on the absolute values of the compact- 
ness parameters, we fix l^ to a constant value in order to im- 
prove the fitting strategy. For l^ in the range 0.5 < k < 100, 
the spectral shape allows to constrain Ith/k to be always of 
the order of 5. Therefore we set Zs = 15, so that the total 
compactness parameter / = Zth + 's is consistent with Eq. [3] 

For the cutoff spectrum, the pure thermal model delivers 
a good description of the data. The best fit (y^/69 = 0.83) 
yields a compactness ratio of Ith/k - ^■^'^-ooh ^ nioderate 
opacity of Tion = 2.40^QQg and a significant ionized reflection 
component (Q/Itt = 0.44 + 0.04 ; ^ = 330+™). The equilib- 
rium temperature of the electrons is found at 32 keV. 

For the excess spectrum, however, the thermal model is not 
appropriate. The best fit is indeed of poor quality, with a re- 
duced chi-square of ;(f^/75 = 1.40. We note that the degra- 
dation results only from the high-energy channels, which are 
responsible for;i' /v = 65/23 of the total j^^^. Obviously, the 
model fails to reproduce the high-energy excess, supporting 
the fact that the emission above 100 keV is linked to non- 
thermal processes. We conclude that a pure thermal model 
is not adequate to explain the observed spectral behavior. 

4.3. Non-thermal acceleration with external Comptonization 

Now we investigate the opposite case, in other words we as- 
sume that the energy injection is purely non-thermal (/nth ^ 0, 
'th — 0). The rest of the model remains unchanged, i.e. we 
consider an external source of soft photons (of fixed compact- 
ness 4 = 15) and we neglect the magnetic field (/^ - 0). Since 
the electrons can partly thermalize through Coulomb colli- 
sions, the equilibrium distribution is hybrid (i.e. thermal/non- 



thermal), even if the whole power is supplied via non-thermal 
injection. The acceleration processes are phenomenologi- 
cally described by the model parameters /nth, Tinj, Tmin and 
Tmax- Here again, we define two sub-models in oder to re- 
duce the number of free parameters. In the first model, we set 
(Tmin, Tmax) = (1-3, 1000), and analyze the spectra in terms of 
/nth and Finj. This configuration is frequently used in the litera- 
ture and will provide a familiar context to discuss our results. 
In the second model, we adopt a more novel approach and 
investigate the eff'ects of a varying maximum energy of the 
accelerated electrons. We keep y^in = 1.3 but fix the spec- 
tral index to the fiducial value Finj = 2.5, while /nth and yn,ax 
are free to vary. For convenience, we call these models the 
ECMl and the ECM2 (for External Comptonization Models), 
respectively. The associated model spectra are computed with 
the Comptonization code eqpair fcoppi,1999,) and the fitting 
results are summarized in Table 3. 

4.3.1. The injection index 

For the cutoff spectrum, the ECMl provides a good fit 
to the data 0^-/68 = 0.84, cf. Figure |5] left). The ob- 
served cutoff shape is well reproduced by a very soft injected 
electron distribution. Our best fit yields a spectral index of 
Finj - 3.95;^Q^Q, which was fixed to this value to determine 
the error on the other free parameters. 

We find a compactness ratio of /nth/4 - 4.35 + 0.08, imply- 
ing a non-thermal compactness of /nth = 65.0 + 1 .0. The opti- 
cal depth of ionization electrons is fitted at Tion = 2.97 + 0.07. 
Because of the soft injected electron distribution, there is only 
very little pair production, increasing the total optical depth 
to tt - 2.99 + 0.07. The observed spectrum strongly re- 
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Figure 5. Joint SPI and PCA spectra fitted with the ECM 1 . The total model is shown by the solid fine, the dot-dashed fine marks the iron fluorescence emission 
from the disc, the dashed line shows the reflection component and the dotted line represents the sum of the Comptonized and reflected emission. The error bars 
are given at the lo" level. 



quires a moderate amount of Compton reflection, with a fitted 
amphtude of Q.j2n - 0.26^qqj and an ionization factor of 
^ = 710 + 120. Freezing Q./2jt to zero leads to a dramatically 
worse fit (F-test probability < 10"^°). 

For the excess spectrum, the best fit (x^/74- = 1.10, cf. Fig- 
ure [5] right) requires a substantially harder power law injec- 

With respect 



tion, with a fitted index of F; 



inj 



2 55 +"-^^5 



25+""- 



to the cutoff spectrum, the compactness ratio increased by 
~15 per cent to /nth/4 - 5-06 +qq". The non-thermal com- 
pactness now yields /nth = 75. 9+ [j and the total optical depth 
increased to tj = 3.87+°°^, which is mainly due to the en- 
hanced production of pairs. The fitted reflection amplitude 
and disc ionization, however, remain stable (0/27r 
; ^ = 600 + 240). 

In order to explore the dependence on compactness, we 
abandoned our fiducial hypothesis /s = 15 and fitted the ex- 
cess spectrum with a variable soft compactness. We find 90 
per cent confidence intervals of 0.3 < U < 250 and 4.5 < 
/nth//s < 5.9, which confirms that the compactness ratio does 
not depend on the total compactness of the source. Although 
our analysis is unable to uniquely determine the total com- 
pactness, the upper and lower limits o n the allowed range ar e 
well constrained by the spectral shape JGierlihski et al .11 9991) . 
The upper bound represents the limit at which, despite an ex- 
tremely soft injected electron distribution (Finj - 4.50+""?), 



the plasma is completely pair dominated (Tit 



10 



-4. 



Tt 



3.0; kTf. ~ 11.5keV). Above this threshold, the growing 
amount of pairs can no longer be balanced by a softer injec- 
tion spectrum, leading to an inaccurate reproduction of the 
thermal peak in the photon spectrum (cf. paragraph 5.2). The 
lower bound corresponds to the limit at which the cooling of 
the high-energy particles is dominated by Coulomb instead of 
Compton losses, leading to an underprediction of the intensity 
of the high-energy tail. A change from /j = 0.3 to l^ - 0.25 



implies a degradation of Ax^ = +5 at 74 d.o.f., showing that 
the threshold is well established. Hence, in the framework 
of an external Comptonization model involving non-thermal 
electron acceleration up to 7,nax - 1000, the total compactness 
of the hard X-ray source can be conservatively constrained by 
means of pure spectral analysis to be 2 < / < 1500. 



4.3.2. The maximum particle energy 

Now, in the ECM2, we freeze Finj to 2.5 but allow for vari- 
ations of the maximum Lorentz factor of the accelerated elec- 
trons. For the cutoff spectrum, we obtain again a very good fit 
to the data ix^ 16% - 0.85). The effects of a much harder slope 
are effectively balanced by a very low cutoff energy of the in- 
jected electron distribution. Indeed, we find yniax = 4.0+,' 3, 
which means that the maximum difference of the kinetic en- 
ergy between the accelerated paricles is about a factor of 10. 
The other parameters are not much affected, namely we find 



/n,h//s = 4.25![;:\'« and T,on - Tt = 3.06 



+0.09 
-0.06- 



Similarly, 



the reflection parameters remain stable {D./27T = 0.23+q"^ ; 

^ = 850!P. 

For the excess spectrum, the fixed injection index is consis- 
tent with the best fit value obtained with the ECMl . Neverthe- 
less, it turns out that there are some differences between the 
results obtained with the two models. Mainly, we notice that 
particle acceleration up to y,^^^ = 21.6+i^^g" is sufficient to re- 
produce the observed high-energy tail. Moreover, a truncated 
electron distribution allows to improve the quality of the best 
fit to ;if^/74 = 1 .00. The total opacity of the plasma is found 
equal to the value obtained with a non-truncated distribution, 
but the pair yield is reduced due to the lack of very energetic 
particles. With respect to the ECMl, the other parameters are 
only marginally affected and remain consistent within the 90 
per cent confidence errors. 
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Table 3 

Best fit parameters of the joint SPI and PCA spectra for the non-thermal models. Fixed paramete rs are indicate d by a star (*) next to their value. The 
temperature of the thermal component is estimated at the equilibrium using Equation 2.8 of lCoppil 09921) and is hence not a fit parameter 



4.4. Non-thermal acceleration with internal Comptonization 

Finally, we consider models in which the observed X/y-ray 
spectra are produced in a magnetized plasma. To emphasize 
the effects of the magnetic field, we study the case where all 
the seed photons are internally generated from synchrotron 
emission (/s - Q, Ib + 0). As in the previously analyzed 
non-magnetic models, we adopt the same different configura- 
tions to phenomenologically describe the acceleration mech- 
anism. For convenience, we call these models the ICMl and 
the ICM2 (for Internal Comptonization Models), respectively. 
To calculate the model spe ctra, we used the new versatile 
Comptonization code belm jBelmont. Malzac. & MarcowithI 
120081) . One of the main differences with eqpair is that belm 
accurately accounts for self-absorbed cyclo-synchrotron radi- 
ation from the sub-relativistic to the ultra-relativistic regime. 
We compared the spectra obtained by both codes for the best 
fit parameters of the non-magnetized models. The relative dif- 
ferences are smaller than 3 per cent at all energy. The fitting 
results are summarized in Table 3. 

4.4.1. The injection index 

First, we employ the commonly used configuration, i.e. we 
fix (ymiiuTmax) = (1.3,1000). A qualitative analysis shows 
that the model spectrum below 30keV is rather insensitive 
to the individual values of /nth, Ib and Finj, but strongly de- 
pends on a combination of all three parameters (see also 
MB09). The situation is similar as in the non-magnetic case 
(cf. section 4.2), although the dependence is slightly more 
complicated. On the other hand, the high-energy spectrum 
(>100keV) is mostly determined by the injection index Finj. 
As a consequence, we use Eq. [3] to fix /nth = 100, while the 
broad band spectra allow to disentangle the degeneracy be- 
tween Ib and Finj ■ 

Due to the detailed treatment of the microphysics in the 
BELM code, real-time fits in xspec are very time consuming. 
In order to make the fitting process more efficient, we tab- 
ulated the model spectrfl The resulting table file has three 
dimensions (corresponding to the parameters Ib, Fi„j and rion) 
and the fits are performed through interpolation between the 
tabulated spectra. 

In order to account for Compton reflection from a 
cold disc, the belm model is convolved with an angle- 
dependent reflection routine based on the pexriv model by 

"1 We used the wftbmd routine (pubhcly available at the heasarc website) 
to create the appropriate FITS file required by the atable model in xspec 



iMagdziarz & Zdziarskil ([l995), but taking into account the 
distortions due to general relativistic effects. The free param- 
eters of the ICMl are hence Ib, Finj, T"ion, Q/27r, ^, the iron line 
energy and the normalization factor between PCA and SPI. 

For the cutoff spectrum, the ICMl provides a good fit to 
the data 01^-/68 = 0.96). In comparison with the ECMl, the 
high-energy rollover is better reproduced. However, as can 
be seen from the residuals in Figure |6] (left), the fit is slightly 
less accurate in the iron line region around 6.4 keV. There are 
also some positive residuals below 5 keV, possibly hinting the 
need for a soft disc component. Since we focus on the spectral 
behavior at high energies, we did not investigate these issues 
any futher. The best fit is obtained with an electron spectral 
index of Finj = 3.47+^^^ The magnetic compactness is fit- 
ted at Ib - 740;'^j4Q, which corresponds to a magnetic field 
strength of B = 2.0+Qg x 10^ G. The plasma is found to be 
of moderate optical thickness, with fitted Tion - 2.26 + 0.18. 
Due to the fast synchrotron cooling of the small number of 
high-energy particles, the produced pair yield is negligible. 
As with the ECMl, the data strongly require ionized Comp- 
ton reflection, with a fitted amplitude and ionization factor of 
D./27T = 0.43 + 0.05 and ^ = 300+J^[J, respectively. 

For the excess spectrum, the ICMl allows again a good de- 
scription of the data (cf Figure |6] right). The best fit yields 
;^'^/74 = 1.06 and the high-energy data constrain the injection 
index to Finj = 2.50;'^qq|. This result is roughly equal to the 
value obtained with the ECMl. However, at equal injected 
electron distribution, the high-energy tail of the ICMl spec- 
trum is slightly steeper. Indeed, contrary to a non-magnetized 
model, the Compton losses compete with the synchrotron 
losses and a significant fraction of the energy radiated by the 
non-thermal leptons is emitted in the optical/UV range rather 
than in hard X-rays. Accounting for the the high-energy tail 
therefore requires more high-energy particles, i.e. a smaller 
Finj. We find a magnetic compactness of /^ = 25.2+3.6, which 
given our assumptions translates to B = 3.76 + 0.14 x 10^ G. 
This means that in the framework of this model, the magnetic 
field strength would drop by a factor of 5.5 during a change 
from the cutoff to the excess spectrum. 

The Thomson optical depth of ionization electrons is fitted 



at Tion = 2.14 



+0.14 



Contrary to the results obtained with the 



ECMl, the Tion values for both spectra are compatible within 
the 90 per cent confidence errors. The reflection routine yields 
Q/27r = 0.44+°°^ and ^ = 260!^,^°, which means that these 
characteristics did not change either. 



The high-energy emission of GX 339^ 



10.0 



Cutoff spectrum 




10.0 



Excess spectrum 



E 
o 
> 10 



0.1 

5.0 
to 
S 0.0 

CO 
Q) 

cc 
-5.0 




V/«¥^V^'" VA i />t i t 



10 



Energy (keV) 



Energy (keV) 



100 



Figure 6. Joint SPI and PCA spectra fitted with the ICMl. The solid line is the total model, the dot-dashed line represents the SSC component while the dashed 
and dotted lines show the reflection component and the iron line emission, respectively. The error bars are given at the Icr level. 



4.4.2. The maximum particle energy 

In the ICM2, we keep /nth = 100 and -/niin = 1-3, but allow 
now for variations of ymax- Instead, we set the injection in- 
dex to Finj - 2.5 and generate a new fitting table which has 
again three dimensions, corresponding to the free parameters 
Ib, Tmax and rion- 

For the cutoff spectrum, we obtain a good fit to the data 
(X^ Id^ - 0.96; shown in Figure [7] left), qualitatively equiva- 
lent to the best fit obtained with the ICMl. Again, this shows 
that the effects of a very soft injection slope can be mimicked 
by a much harder injected distribution, but truncated at a cer- 
tain particle energy. We find yn,ax = 15.2^j-^, which is sig- 
nificantly higher than the fitted maximum particle energy in 
the ECM2. The inferred magnetic compactness is reduced 
with respect to the ICMl, that is to say h - 280+3^, which 
for a medium of typical size /? = 30 Rq corresponds to B = 
1.25^Q4 X 10' G. The other parameters are not much affected, 
namely we find Tion - tj 
component we infer Q./2n 

For the excess spectrum, we note that the injection index 
is fixed to the best fit value obtained with the ICMl. How- 
ever, allowing for ymax to vary (i.e allowing for a truncated 
electron distribution), the fit may be improved considerably. 
For ynjax = 190^[q°, we obtain the best description of the ex- 
cess spectrum, with a reduced jt'^ of;^'^/74 = 0.92 (cf. Figure 
l2]right). Compared to the non-magnetized model, the maxi- 
mum electron energy is again found to be much higher This is 
expected since in the magnetic models, the seed photons have 
a lower average energy than in the non-magnetized models. 
Indeed, in the ICMs, the high-energy photons are produced 
from single Compton scattering off the synchrotron emission 
(£5 ~ 0.01 keV) while in ECMs the seed photons originate 
from the disc emission (£5 ~ 1 keV). Therefore, in order to 
upscatter the seed photons to 200 ke V, the electrons must have 



12 in the ECMs and y ~ 120 



2.30+°|^ and from the reflected 

0.43!0:04and^ = 300!j20. 



averaged Lorentz factors of y 
in the ICMs. 

The magnetic compactness remains equal to the value in- 
ferred with the ICMl, namely Ib - 27.3^34. A transition 
from the cutoff to the excess spectrum thus requires a factor 
3.2 decrease of the magnetic field strength. The opacity from 
ionization electrons yields Tion = 2.49 + 0.18, while the total 
optical depth is found to be tj = 2.73+0.18. This shows again 
that the ionization opacity does not change between the two 
spectra. Finally, the fitted reflection amplitude and ionization 
parameter are found to be consistent for both spectra regard- 
less of the acceleration model (cf. Table 3). We thus conclude 
that in the framework of a strongly magnetized medium, the 
constraints to these parameters are very robust. 

5. DISCUSSION 

The SPI observations showed that during the bright 
hard state of the 2007 outburst, the highest energy emis- 
sion (>150keV) of GX 339-4 was variable. While the 
spectral shape at lower energies (4-150 keV) remained 
more or less constant, we detected the significant appear- 
ance/disappearance of a high-energy tail. The strength of this 
hard tail, varying on a time scale of less than 7 hours, is found 
to be positively correlated with the total X-ray luminosity of 
the source and enables interesting constraints to the physi- 
cal processes which could be responsible for the high-energy 
emission. 

5.1. The pure thermal model 

The clear detection of a cutoff energy in the cutoff spectrum 
indicates that the Comptonizing electron distribution is quasi 
Maxwellian. Thus, it is not surprising that the spectrum can 
be explained by assuming only thermal heating of the plasma. 
As suggested by the advection dominated accretion flow mod- 
els, this could be achieved through Coulomb interactions with 
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Figure 7. Joint SPI and PCA spectra fitted with the ICM2. See Figure|6]for the description of the dift'erent model components. The eiTor bars are given at the 
\iT level. 



a thermal distribution of hot protons. The temperature of the 
protons can be estimated from the thermal compactness, the 
electron temperature and the optical depth of the plasma (cf. 
formula (4) in MB09). For /th = 100, we infer a proton tem- 
perature of the order of 1 MeV. This is of the same order of 
magnitude than the proton temperature estimated in MB09 
for the canonical hard state of Cygnus X-1, which in com- 
parison to the hard state of GX 339-4 analyzed here shows 
a hotter electron plasma (kT^ ~ 85 keV) but a lower com- 
pactness (/th ~ 5). As mentioned in MB09, proton tempera- 
tures of the order of 1 MeV are significantly lower than what is 
expected in typical two-temperature accretion flows, namely 
kTi > 10 MeV. 

In any case, pure thermal heating is not enough since it 
is not able to explain the appearance of the observed high- 
energy tail. We conclude that either the hard excess is inde- 
pendent from the thermal component, in which case its origin 
is located outside the innermost regions, or that both compo- 
nents are linked, in which case at least some level of non- 
thermal heating is required. 

On the other hand, both broad band spectra can be suc- 
cessfully explained by models involving only non-thermal 
electron acceleration. This suggests that the Comptonizing 
medium in hard states could be powered by the same non- 
thermal mechanisms that are believed to power the accretion 
disc corona in soft states. Depending on the nature of the 
plasma (magnetized or not) and the origin of the seed pho- 
tons, these issues are discussed in the next paragraphs. 

5.2. The non-magnetic case 

In the framework of a non-magnetized model, the vari- 
ability of the high-energy spectrum can be explained by 
changes in the properties of the involved acceleration pro- 
cesses. Namely, the fits suggest a small variation of the total 
power supplied to the plasma along with a significant varia- 
tion of either the spectral index (in the ECMl) or the cutoff 



energy (in the ECM2) of the non-thermal electron distribu- 
tion. 

In the ECMl, a change from the cutoff to the excess spec- 
trum requires that the spectral index typically drops from 
Finj = 4.0 to 2.5. This implies that the average energy of the 
accelerated particles <Einj> rises from 1.0 to 1.8 MeV In the 
ECM2, at constant spectral index, the best fits show that rel- 
atively low maximum Lorentz factors are sufficient to repro- 
duce the data. In this case, the appearance of the high-energy 
excess requires an increase of ymax from 4 to 22, which im- 
plies that the average energy of the accelerated electrons rises 
accordingly, from 1.0 to 1.5 MeV. In both cases, such vari- 
ations can be explained by the possible non-stationarity of 
the inherent acceleration mechanisms. For instance, in the 
framework of shock acceleration, the propertie s of the ac- 
celerated particles depend on the shock strength dWebb et al.l 
|T984; Spitkovskv 2008). Acceleration by reconnection de- 
pends on several physical pa rameters such as t he local geome- 
try of the reconnection zone (iZenitani and Hoshino 20071) and 
the number of reconnection sites (if the particles are acceler- 
ated stoc hastically by successive acceleration events in difFer- 
ent sites jAnastasiadis et al.|[T997t iDauphin et al.ll2007b ). AH 
these properties may undergo variations with time and could 
hence explain the observed variability. 

Our results are con sistent with those of 
iGierlinski & Zdziarskil (120051) . who studied the energy 
dependent variability of non-magnetized Comptonization 
models in response to varying physical parameters. Although 
the variations of Finj and y^ax could not explain the observed 
rms spectra of XTE J 1650-500 and XTE J 1550-564, they 
found that changes in these parameters produce strong 
variations in the X-ray spectrum above 50 keV, which is what 
is required here to reproduce the present data. 

The spectral analysis suggests that a transition from the cut- 
off to the excess spectrum additionally requires a 15 per cent 
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Figure 8. Photon spectra (left panel) and particle distributions (right panel) in the cutoff and excess states, obtained with the belm code for the best fit parameters 
of the ICM2. The photon spectra include neither absorption nor reflection. T(p) is the Thomson optical depth per unit momentum. For each particle distribution, 
the vertical dotted lines show the energy chosen to define the thermal and the non-thermal populations (cf. section 5.3). 



increase of the total power supplied to the plasma. Consid- 
ering a spherical medium of fixed radius and constant illumi- 
nation from the accretion disc, this increase is consistent with 
the 14 per cent difference in the observed 4-500 keV luminos- 
ity. In addition, the total optical depth of the plasma is found 
to increase by about 30 per cent, which is mainly due to the 
enhanced pair production occuring in the harder acceleration 
regime. 

Motivated by size and luminosity estimates, we assumed a 
constant soft photon compactness of 4 - 15. As mentioned 
earlier, the above results are only weakly dependent on the ex- 
act value of /j. Very high- or very low values of the illumina- 
tion compactness, however, turned out to be inconsistent with 
the observed spectra. Indeed, since the ratio /nth/4 is robustly 
constrained, strong illumination requires an efficient accel- 
eration mechanism (i.e. large /nth) to reproduce the spectral 
slope at lower energies (< 20keV). Consequently, this gener- 
ates very energetic radiation which produces large amounts of 
e e^ pairs through photon-photon annihilation. This, in turn, 
reduces the equilibrium temperature of the plasma since more 
particles have to share the same amount of energy. Hence, 
above a certain soft photon compactness, the equilibrium tem- 
perature will be too low to be consistent with the observed 
thermal peak of the spectrum. 

Reciprocally, since the thermal part of the electron distribu- 
tion is roughly determined by the balance between accelera- 
tion and Compton cooling, decreasing /, at constant /nth//s has 
no significant effect on the lower energy part of the photon 
spectrum. However, a major fraction of the high-energy tail 
results from the Comptonization by mildly relativistic parti- 
cles (2 < y < 10). Below a certain soft photon compactness, 
the cooling of these mildly relativistic particles is no longer 
dominated by the Compton losses but by Coulomb interac- 
tions with the lower-energy thermal electrons. Thus, decreas- 
ing /s at constant /nth//s provides a weaker acceleration rate 
while the cooling remains constant. This reduces the inten- 
sity of the high-energy tail up to the point where the model 
predictions are no longer consistent with the data. 

In conclusion, independently of any geometric argument, 
we obtain conservative bounds to the total compactness of 
the X-ray emitting plasma, i.e. 2 < / < 1500. These limits 
are consistent with the estimates derived from geometric argu- 
ments, but unfortunately not very constraining. Anyhow, the 



robustness of the fitted compactness ratio /nth//,s - 4.5 allows 
to conclude that the luminosity of the cold disc represents at 
most ~ 20 per cent of the luminosity of the Comptonized com- 
ponent, possibly much less if the plasma is magnetized. 

5.3. The magnetic case 

Usi ng the new code belm jBelmont. Malzac. & MarcowithI 
I2008h . we showed that the hard X-ray behavior of GX 339- 
4 in the bright hard state can be explained by assuming 
pure non-thermal electron acceleration and subsequent Comp- 
tonization of the self-consistently produced synchrotron pho- 
tons. The model requires no incident radiation from the ac- 
cretion disc and assumes constant power injection into the 
magnetized plasma. As in the non-magnetic case, the spectral 
variability can be mimicked by two different configurations 
of the acceleration model, involving either a variable power 
law slope (in the ICMl) or a variable maximum energy (in 
the 1CM2) of the injected electron distribution. 

In principle, these models allow to estimate the averaged 
magnetic field strength of the plasma. However, since the fits 
provide precise constraints only for the ratio IbIU the uncer- 
tainties on the total compactness (cf. equation Q) are pro- 
jected to the estimate of Ib- To discuss our results, we express 
Ib as a fraction of the magnetic compactness at equipartition 
with the radiation field /g^. As we have the approximate de- 
pendence /br (^ / X (1 + tt/3) (cf. equation (8) in MB09), the 
ratio IbIIbr does not depend on the uncertainties regarding the 
source size and distance and is therefore a good indicator to 
quantify the importance of the magnetic processes. In addi- 
tion, /s = was assumed to study the physics of a strongly 
magnetized medium, but it can not be excluded that both the 
synchrotron flux and a soft disc component contribute to the 
cooling of the non-thermal particles. If the medium is ad- 
ditionally illuminated by cold disc photons, less synchrotron 
cooling will be required to reproduce the slope of the lower 
energy spectrum, implying that the fitted values of the mag- 
netic compactness are in fact conservative upper limits. In the 
ICMl, we infer IbIIbr < 18 from the cutoff and Ib/Ibr < 0.58 
from the excess spectrum. In the 1CM2, the fitted magnetic 
compactness for the cutoff spectrum is lower, namely we find 
h/h^ ^ 6.0 and Ib/Ibr ^ 0.60, respectively. 

As mentioned above, the magnetic models do not require 
any disc blackbody photons to produce the observed 4- 
500 keV spectra. If the accretion disc extends down close to 
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the black hole (as requested by the accretion disc corona mod- 
els), the paucity of soft disc photons can be explained if the 
Comptonizin g coronal material is outflow ing at a mildly rela- 
tivistic speed (iBeloborodovl 1 9991: iMalza c et al. 2001). Using 
the formulae (5) and (7) of lBeloborodovl (I1999D . we estimate 
that bulk velocities of at least 0.6 c are required for the cooling 
of the corona being dominated by synchrotron self-Compton. 
Comptonization off a dynamical corona may then blueshift 
the emerging spectrum, but these corrections remain moder- 
ate at ~0.6 c and have not been included in the spectral fits. On 
the other hand, the disc may as well be truncated (as requested 
by the hot flow models), since the data do not explicitly re- 
quire relativistic smearing of the reflection features. Thus, as 
long as the particle acceleration is essentially non-thermal, the 
ICMs may apply to both geometries. 

If the electron cooling is dominated by synchrotron self- 
Compton (k/lnth ^ 1/5 from the non-magnetic models), 
the fits suggest that the magnetic field strength is roughly in 
equipartition with the radiation field. From a qualitative fit 
of the canonical hard state spectrum of Cygnus X-1, MB09 
constrained the magnetic energy density to be strictly below 
equipartition (Ib/Ibr < 0.3). As a consequence, even if our 
results are consistent with the results obtained for Cygnus X- 
1, the magnetic field could play a more important role in the 
physics of GX 339^, at least in the cutoff state. 

In both magnetic models, a transition between the two spec- 
tra can be explained by the variations of only two parameters. 
In any case, the magnetic compactness Ig needs to be vari- 
able, along with either the injection slope (in the ICMl) or the 
maximum energy of the accelerated particles (in the 1CM2). 
All other fit parameters are found to remain constant within 
the 90 per cent confidence errors. To reproduce the cutoff to 
excess transition with the ICMl, the magnetic field has to de- 
crease by a factor of 5.4 and the injection slope drops from 
Finj - 3.5 to 2.5, while in the 1CM2, the magnetic field de- 
creases by a factor of 3.2 and the maximum energy increases 
fromy^ax = 15.2 to 187. 

Regardless of the precise acceleration process involved in 
the accretion flow, the inferred change in the magnetic field 
strength is expected to have an impact on the cutoff energy of 
the accelerated electron distribution. Indeed, the maximum 
energy of the particles is achieved when the energy losses 
become larger than the gains. In the magnetized models in- 
vestigated in this paper, the energy losses of the most ener- 
getic particles are dominated by synchrotron cooling (Comp- 
ton and Coulomb cooling are significantly smaller), which ob- 
viously depends on the magnetic field strength. In the frame- 
work of diffusive shock acceleration for instance, it has been 
shown that when synchrotron losses are included, the max- 
imum Lorentz factor of the accelerated relativistic electrons 
satisfies ymax °c l/(KB^), w here K is the diffusion coefficient 
(see e.g. I Webb et al.l [19841 or|Marcowith & Kirk 1999). If K 
is constant, this predicts a maximum energy of yniax °^ l/B^. 
However, depending on the assumptions made to describe the 
acceleration mechanism, the diffusion coefficient can depend 
both on the particle energy and the magnetic field strength: 
K{y, B). In the frequently used Bohm limit, it is assumed that 
K oc y/B, implying that the maximum Lorentz factor of the 
accelerated particles is expected to follow y„,ax '^ B"'''^. Al- 
though our results are not consistent with the Bohm predic- 
tions, the overall behaviour remains that the cutoff energy de- 
creases with an increasing magnetic field strength. Moreover, 
other acceleration mechanisms, such as magnetic reconnec- 
tion for instance, may give different predictions. 



As mentioned in the previous section, the injection index 
is not universal and may undergo variations in response to 
changing physical conditions in the acceleration region. How- 
ever, while mere synchrotron cooling implies that y^ax and 
B are necessarily anti-correlated, the physical connection be- 
tween the variations of Finj and the magnetic field strength 
remains less obvious. Hence, in the context of a magnetized 
plasma, the model with variable y^n^ is slightly favored. Any- 
how, the variability of the acceleration mechanism was mod- 
eled by the variations of only one parameter (either the in- 
jection index or the cutoff energy), although it cannot be ex- 
cluded that both parameters undergo simultaneous variations. 
We emphasize however that such models are more compli- 
cated for they imply more varying parameters and would not 
significantly improve the quality of the fits. 

In the context of the 1CM2, we investigated the impact of 
the parameter variations on the electron distribution and the 
involved radiation mechanisms. Although the lower energy 
part of the model photon spectra is similar in shape, the under- 
lying electron distributions are quite different (cf. Figure [8]). 
Indeed, the inferred variations of B and ymax not only change 
the non-thermal part of the distribution but also the temper- 
ature of the thermalized component. Figure |9] compares the 
contributions from the various radiation mechanisms to the 
model photon spectra, showing separately the Comptoniza- 
tion off thermal and non-thermal particlefl As expected, the 
parameter variations strongly increase the ratio of the non- 
thermal to the thermal component, resulting in the production 
of the high-energy tail. 

In conclusion, the 1CM2 provides a framework for a sim- 
ple, physically motivated interpretation of the data, showing 
that the spectral variability could be triggered by the varia- 
tion of only one single parameter, namely the magnetic field 
strength. Assuming a spherical medium of radius R - 30 Rg, 
a black hole mass of M = 13 M© and a distance of d - d, kpc, 
we infer a factor 3 variation of the magnetic field strength, be- 
tween B < 3.9+}-^ X lO'-G and B < 1.25+°-^^' x lO^G. If the 
fitted average spectra provide a good estimate of the individ- 
ual spectra from the single science windows, the time scale 
of this evolution is at m ost of the order of hours. Using the 
standard o'-prescription dShakura & Sunvaevlll973h . the vis- 
cous time scale of the accretion disc at a radius R is given 
by fvisc = « ' (HIRY" Ik, where H/R is the aspect ratio of 
the disc and Jk the Keplerian period. Using the lower limit 
a > 0.01 (quiescent disc) and H/R ^ 0.1 (thin disc), we 
obtain fyisc < lOmin for the typical source size R = 30 Rq. 
The viscous timescale of ADAF-like models is much shorter. 
Therefore, even if the geometry of the X-ray emitting region 
and its dynamical evolution remain uncertain, global changes 
on time scales of hours are not unrealistic. 

6. SUMMARY AND CONCLUSION 

We presented an analysis of the high-energy emission of 
GX 339-4 in a luminous hard state. With respect to the stan- 
dard cutoff shape of the hard state spectrum, the 25-500 ke V 
INTEGRAL/SPl data revealed the appearance of a variable 
high-energy excess. The intensity of this hard excess seems 
to be positively correlated with the total X-ray luminosity and 
the associated time scale is shorter than 7 hours. We explored 
the possible physical origins of this variability through an ex- 

^ The low energy part of the particle distributions were fitted with a 
Maxwellian and particles of energy one order of magnitude larger than the 
infen'ed temperature were considered to be non-thermal. 
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respectively. 



tensive analysis of two averaged spectra, one showing the typ- 
ical cutoff shape and one showing this prominent high-energy 
excess. 

We used simultaneous RXTEfPCA data to extend the spec- 
tral coverage down to 4 keV and fitted the broad band spectra 
with a variety of physical Comptonization models. Models 
involving only thermal heating can be ruled out since they 
are not able to reproduce the high-energy tail. This feature 
thus confirms that in luminous hard states, the Comptoniz- 
ing plasma of GX 339-4 contains a fraction of non-thermal 
particles. Models involving only non-thermal electron accel- 
eration, on the other hand, showed that the thermal part of the 
spectrum can be produced by an initially non-thermal (power 
law) distribution, which rapidly thermalizes under the effects 
of synchrotron self-absorption and/or e-e Coulomb collisions. 

The relatively good signal to noise ratio of the high- 
energy channels (>150keV) allowed to derive meaningful 
constraints to the model parameters. Depending on the na- 
ture of the plasma (magnetized or not), the transition between 
the two averaged spectra requires the variations of at least two 
parameters. We found that a magnetized medium subject to 
non-thermal electron acceleration provides the framework for 
a simple and physically consistent interpretation of the data. 
Indeed, we showed that the spectral variability can be trig- 
gered by the only variations of the magnetic field, implying a 
subsequent variation of the maximum energy of the acceler- 
ated particles. 

The quantitative constraints derived from this model yield 
a very conservative upper limit on the average magnetic field 
strength in the Comptonizing plasma and suggest that in the 
bright hard state, the magnetic energy density could reach 
equipartition with the radiative energy density. 

In conclusion, the presented results suggest that magnetic 
processes are likely to play a crucial role in the production 
of the high-energy emission of GX 339-4. In luminous hard 
states, the Comptonized emission could originate from a mag- 



netized corona essentially powered though non-thermal parti- 
cle acceleration, similarity to what is believed to happen in 
soft states. The hard X-ray emission in both spectral states 
may therefore be the consequence of a common physical phe- 
nomenon. 
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